Molecular switch controlled by pulsed bias 

voltages 

Velimir Meded, Alexei Bagrets, Andreas Arnold, and Ferdinand Evers 

June 25, 2009 



*PD Dr. Ferdinand Evers 
Institut fiir Nanotechnologie 
Forschungszentrum Karlsruhe 
D-76021 Karlsruhe, Germany 
Institut fiir Theorie der Kondensierten Materie 
Universitat Karlsruhe 
D-76128 Karlsruhe, Germany 
E-mail : Ferdinand .Evers@int.fzk.de 
Dr. Velimir Meded, Dr. Alexei Bagrets 
Institut fiir Nanotechnologie 
Forschungszentrum Karlsruhe 
D-76021 Karlsruhe, Germany 
Dr. Andreas Arnold 

Institut fiir Theorie der Kondensierten Materie 
Universitat Karlsruhe 
D-76128 Karlsruhe, Germany 



1 



Keywords: Density functional calculations, Molecular electronics, Molec- 
ular switches. Molecular wires. Molecular mechanics 

Abstract 

It was observed in recent experiments that the current-vohage 
characteristics (IV) of BPDN-DT (bipyridyl-dinitro oligophenylene- 
ethynylene dithiol) can be switched in a very controlled manner be- 
tween "on" and "off" traces by applying a pulse in a bias voltage, 
^ias- We have calculated the polaron formation energies to check a 
frequently held belief, namely, that the polaron formation can explain 
the observed bistability. Our resuhs are not consistent with such a 
mechanism. Instead, we propose a conformational reorientation. The 
molecule carries an intrinsic dipolc moment which couples to T4iias- 
Ramping V^as exerts a force on the dipole that can reorient ("rotate") 
the molecule from the ground state ( "off" ) into a metastable configura- 
tion ("on") and back. By elaborated electronic structure calculations, 
we identify a specific path for this rotation through the molecule's 
conformational phase space. We show that this path has sufficiently 
high barriers to inhibit thermal instability but still the molecule can 
be switched in the voltage range of the junction stability. The theo- 
retical IVs reproduce qualitatively the key experimental observations. 
We propose, how the alternative mechanism of conductance switching 
can be experimentally verified. 
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1 Introduction 



During the past few years, experiments in the field of "Molecular Electronics" 
have experienced significant improvements. This made possible the observa- 
tion of Coulomb blockade [Ij and Kondo effects |2] in molecular transistor 
geometry. Improved statistical methods to access a conductance of molec- 
ular junctions have been developed [31 H]. A spectacular manifestation of 
molecular low energy excitations has been detected [5] in inelastic current 
spectroscopy. 

Aiming for potential applications, Molecular Electronics holds the vision 
that functional devices like memory elements, switches, transistors etc. may 
be realizable by designing suitable molecular complexes. In this spirit, a 
considerable effort to design and test molecular systems with a controlled 
switching behavior has been made P, [7] . The challenge here is a ubiquitous 
one in the field of Nanotechnology, namely how to build and manipulate a 
device or a material with full control down to the atomic level. Therefore, 
one may expect that a solution in one subfield will cross-fertilize exciting 
developments in many others. 

The proposed prototypes for single molecule memory elements can be 
sorted according to which physical degree of freedom underlies the switching 
bistability. The three principle categories are: (I) switches with different 
(meta-)stable charge configurations [HI El [IHl [HI [121 [13] ; (H) conformational 
switches with molecules exhibiting different stable isomers [T^ [T5l [T6l ITf] ; 
(III) conformational switches with molecular reorientations against the con- 
tact atoms [L8\ . Once it is known how to operate a single molecule, arrays of 
molecules can be addressed as well, which eventually may form an entirely 
new class of molecular hybrid-materials. 

It is very encouraging, that in recent experiments switching has been 
successfully demonstrated by Blum et al. [I9], Keane et al. [20] and Lortscher 
and Riel |2Ij. They have been using a molecular wire of the "Tour-type", 
BPDN-DT, see Fig. [H which was investigated in several experimental and 
theoretical works, before [ISl [221 [23] • This particular set of experiments 
has two intriguing features: first, switching is established in a two-terminal 
device, just employing a pulse in the bias voltage, and second the effect is 
very stable, so that it has even been observed at room temperature. 

Specifically, the BPDN-DT wire shows the following characteristic behav- 
ior in experiment [21]: with a very slow (adiabatic) change of the bias, Vbias, 
the current voltage characteristic (IV) is fully reversible for up-sweeping and 
down-sweeping. However, with faster sweeping rate the adiabatic regime is 
left and a hysteretic behavior is seen in the IV if Vbias is tuned beyond ~ IV. 
The new IV-curve signahzes the existence of a second (meta-) stable state of 
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the molecular junction. One returns to the original curve only by repeating 
the procedure with inverted bias. 

The physical origin of bistability in these experiments has not yet been 
identified up to now. He et al. [23] attribute it for their experiment in elec- 
trochemical environment to a change in the oxidation state of the molecule. 
A polaron caused bistability was put forward by Galperin et al. [21] also as a 
possible explanation of switching observed in the break junction experiment 
by Lortscher et al. [21] who investigated a freely suspended molecular bridge. 
On the other hand, Keane et al. [20] argue that the formation of a polaron 
is not a likely occurrence in their electromigrated break junction experiment. 
They suspect that rather a bias driven modification of the contacts may be 
responsible. 

Identifying a possible switching mechanism operating for Tour-type mole- 
cules is the aim of this work. Our calculations suggest the following picture: 
the BPDN-DT molecule in vacuum or in between two contacts, Fig. [2], is not 
susceptible to charging. Unless it is stabilized by counter-charges that are 
located very close to the bipyridine unit, any excess charge recombines with 
its images on the electrodes (Fig. E]). This makes polaron-formation (class I) 
unlikely in the experiments [HI [20[ [21] in support of an earlier claim [20j . 

However as we shall demonstrate, the key property for switching, bistabil- 
ity, can originate from a rotational degree of freedom associated with BPDN- 
DT molecule contacted to electrodes. The inertia that is driving the rotation 
results from the action of the electrical field (associated with Vbias) on the 
dipole moment of the dressed bipyridine unit. Two stable configurations can 
be reached by a double axis vr-rotation, see Fig. H] (classes II & HI). Our 
calculations suggest, that the rotation can be performed in such a way, that 
the rotational barrier is high enough — so both states are stable against tem- 
perature — but at the same time this barrier is low enough so that one state 
rotates into the other under bias voltages ~ IV (Fig. [5]). The theoretical 
IV curves (Fig. [6]), which one obtains for either state, closely resemble the 
experimental findings [IHl [201 121] ■ 

2 Model and Method 

The model setup used for the DFT transport calculations is depicted in 
Fig. [21 The Aui4-clusters mimic the electrodes. Computations based on the 
density functional theory (DFT) have been performed with the real-space 
based package TURBOMOLE [53]. Optimized basis sets of triple-C quality 
including polarization functions have been used (exception: IV with a split- 
valence basis set of double-^ quality) [211 [2Z] • The exchange correlation (XC) 
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functional BP86 [2H1 [22] was employed. Charging analysis and relaxations 
were checked against a hybrid functional, B3LYP {3U] . 

3 Absence of polarons 

In order to investigate the possibility of polaron formation, we have consid- 
ered the core region of the molecule, only, replacing the Aui4 cluster by an 
H-atom. We compared the electronic structure for the uncharged (Q=0) and 
the charged species {Q= ± 1 in units of the electron charge |e|). Energies for 
the highest occupied (lowest unoccupied) molecular level, HOMO (LUMO), 
are given in Fig. El We first discuss the molecule in the gas phase, "free" 
molecule. The data indicate, that the LUMO of the positively charged ion 
{Q= + 1) has an energy -Elumo(Q= + 1) ~ — S.leV. Upon approaching a 
metal surface, -Elumo(<5= + 1) will increase due to the interaction with the 
image charge. Assuming, that the excess charge distribution is arranged 
about the center of the molecule with an extension much smaller than the 
molecule length, one can give an estimate of the energy change by employing 
an image charge analysis. We obtained an energy -E'lumo(Q= + 1) ~ — 7.2eV 
for the LUMO after screening. Since it is certain [31], that the workfunction 
of any uncharged Au-surface is above — 6eV and below — 5eV, one can con- 
clude that the positively charged ion will be neutralized as soon as it makes 
contact to a Au-surface. By a completely analogous argument for the HOMO 
of the anion (see Fig. [3]), we deduce that also the molecule with Q= — 1 is 
unstable near a Au-surface. 

All calculations have been performed with a fully relaxed molecular struc- 
ture; the charged molecule undergoes a conformational change, which mostly 
affects the N02-groups. However, our calculations do not give any indica- 
tion, that a geometrical deformation could stabilize the ion in the vicinity of 
the electrodes, see Fig. [3l More precisely, the energy change of the anion's 
HOMO (cation's LUMO) upon structure relaxation is not sufficient in order 
to induce a level crossing with the Fermi energy. This implies, that the con- 
tacted BPDN-molecule does not exhibit a polaron inside a trivial vacuum 
bounded by the two electrode plates. This is the case in the break junction 
experiments [20], [21]. Our conclusion does not necessarily apply to experi- 
ments in an electrochemical environment [23] if this has a sufficiently large 
dielectric response. 
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4 Rotations and bistability 



After we have ruled out charge degrees of freedom as a hkely origin for bista- 
bihty, we now turn our attention to mechanisms related to charge neutral 
conformational changes. Since the control molecule, Fig. [T|, did not show 
hysteresis [T^ [21], bistability should involve a (charge neutral) conforma- 
tional modification of the BPDN-functional unit. An important second re- 
quirement is that an external electric field (realized via Hias) must be able to 
address the putative degree of freedom. Since the monopole moment of the 
functional unit remains zero, it should be the force on its dipole moment, po, 
generated by the N02-groups that pushes one configuration into the other 
and back; according to our calculations (employing BP functional) po ~ 3.0 
Debye. Unlike the BPDN-DT, the reference molecule (BP-DT) enjoys the 
inversion symmetry and, therefore, has a vanishing intrinsic dipole moment. 
We take this an explanation why its IV does not exhibit hysteretic behavior. 

By this reasoning, one is led to look for rotations that take either the 
entire molecule or at least its BPDN-unit over from a stable ground state 
configuration (GC, "off") to another (meta-)stable one (MC, "on"). Any 
such rotation corresponds to a path in the atomic configuration space of the 
(extended) molecule. Apart from bistability, the optimal path has to satisfy 
a number of constraints in order to be a viable candidate for realizing a 
molecular switch, (i) The energy barrier between GC and MC, AE*, should 
exceed temperature, AE* ^ T, to avoid uncontrolled thermal switching, (ii) 
In external fields the optimal path should have a continuous deformation so 
that GC at forward bias, say, and MC at reverse bias become unstable in 
the sense that AE* < T. The instability should be reachable at "switching" 
voltages. Von, still tolerable for the molecular junction, in experiments Von ~ 
1 Volt, (iii) GC and MC should have two IV curves, Jon and Iqs, which can 
be discriminated from each other. 

Giving typical switching voltages Km ~ 1 Volt, one can estimate the 
energy barrier AE* separating two states at zero bias. Assuming that the 
molecule (length d^25 A) forms an angle ca. 45° with the bias electric field, 
we have £on—V^Von/d. Once a dipole, po, is flipped in the bias field, the 
system's energy gain is Ai^dipoie— 2po^^on/ V^, which should be of order AE* 
if energy mismatch between GC and MC is small compared to AE*. Using 
Po — 3 Debye, one arrives to an estimate AE* ~ Aii^dipoie~50 meV. 

After an extensive search, we have found a family of paths which is best 
described as a consecutive one parameter (angle 6) rotation about two axes, 
see Fig. HI The family parameter is the (average) angle, that is formed 
between the longitudinal axis of the molecular wire and the (hypothetical) 
surface normal of the electrode. 
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Consider a variation of ground state energy of the system, AEtot{0), with 
.rotation angle 6. As may be inferred from Fig. (trace OV), the example path 
with ^^70° satisfies condition (i), since the barrier between GC (^=0) and 
MC {e^ir) is AE*^AEtot(vr/2)~100meV^1160K which is in rough agree- 
ment with the above estimate. (Notice, that the double axis rotation is 
crucial to obtain consistent energies. Namely, unbalanced single rotations, 
e.g. about the S-S axis or the C-C axis near the functional unit, give much 
higher barriers, ~ 200 meV, and thus are ruled out.) 

Furthermore, the evolution of the traces depicted in Fig. [5] under the 
applied homogeneous electric field, Sh, shows that for a path with ^^70° 
also condition (ii) is met: the value of at which either GC or MC becomes 
unstable corresponds to an effective voltage iKnl^S -^ 4 Volt. 

Finally, we can show that also (iii) is fulfilled. To this end we have 
determined the IV curve for the situations GC and MC, again at $?a70°, with 
our home made transport package employing the non-equilibrium Green's 
function method [32] . 

5 IV-characteristics 

The theoretical IV-curves, Fig. [6l clearly displays two different IV traces 
which allow to read out the state of the molecular bit, "on" or "off". They 
exhibit key features, which now we discuss: 

(a) A step-like increase in the current is observed (near IV). It is mainly 
due to the molecular orbital HOMO* entering the voltage window. 

(b) Under the transformation 

X: Hias^-Hias,/^-/ (1) 

both IVs are nearly invariant, see insets Fig. [HJ Eq. ([T]) implies that dl/dV 
is an even function of Vbias- Generally speaking, X— invariance is an exact 
property of non-interacting (i.e. non- polarizable) electron systemsQ In inter- 
acting systems, charge localized on the molecule is redistributed as a response 
to Vbias- The polarization produces a change in the effective potential that 
feeds back into the current carrying orbitals thus giving rise to second order 
effects, Vj^^ias, on the current. Details of the redistribution depend strongly 
on contacts and the orientation of the molecule in the bias induced E-field 
(in particular its sign). In the present case the molecular junction has an 
approximate inversion symmetry (see Fig. HI) and therefore violation of X is 
weak. 

^In our discussion we ignore any structure in the density of states, g{E), of the elec- 
trodes; i. e. we assume that g{E) is independent of energy. 
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(c) It is important for understanding the theoretical IV to notice, that 
symmetry violation is stronger with "on-" than with the "off-" trace (inset 
Fig. [6]), so that Jon > /off at positive Vbias but |/on| < |/off| at sufficiently large 
reverse biases | Vbias | > K (dashed trace in Fig. [6] is always above the black 
one). For the ^-rotations, Fig. HJ the difference relates to a slightly modified 
contact geometry. 

(d) The contact difference between "on-" and "off-" states also manifests 
itself in the zero bias conductance; according to our calculations Goff=0.0016 
and Gon=0. 00175, so Goff ^ Gon- Together with (c), the inequality implies 
that there must be an intersection of the IVs at a voltage V^, which is a detail 
observed in our data. Fig. [6l upper inset. 

6 ^-invariance of <l>-junctions and IV 

For a more general class of molecular $— junctions, which do not share the 
(approximate) symmetries of the extended molecule. Fig. |U X— invariance 
can be strongly broken if polarization effects become large. In this case the 
hysteretic features of the IV will be even more pronounced. Indeed, there is 
an intrinsic residual symmetry of any ^-junction, 

/on(Vbias) = —/off( — Vbias), (2) 

which in reality will be violated but only due to contacts. If the impact of 
contacts is not too strong, then the approximate validity of Eq. ([2]) ensures, 
that the junction's IV-curves closely resemble the switching characteristic of 
Fig. [HI (main plot) in the sense that there is a well defined upper ("on") and 
a lower ("off") curve. 

7 Discussion 

The theoretical IV, Fig. [H] shares all essential qualitative features with the 
experimental result [21]. For three reasons, a detailed quantitative compar- 
ison of absolute values for currents and the characteristic voltages, Vc and 
Von, appears quite difficult. (1) Our theoretical analysis implies that the 
angle $ determines the switching voltage. Von-, since it controls the effective 
force (~ sin $) that acts on the intrinsic BPDN-dipole moment facilitat- 
ing the molecular rotation. Since details of the atomic conformation of the 
molecular junction realized in experiments are not known, $ remains largely 
unspecified. In fact, our atomistic modeling of the contact. Fig. [21 can be 
only a crude caricature of experimental reality and therefore parameters like 
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$ are not necessarily very well defined; in addition their statistics may fluc- 
tuate between different experiments [121 ED ED]- (2) Another uncertainty 
relates to the fact that in experiments the molecule may not be in a fully 
relaxed conformation. Our preliminary calculations show that depending on 
whether or not molecular geometries are relaxed, the ii^tot(6')-traces depicted 
in Fig. O vary, if only in a quantitative way. Moreover, our switching path. 
Fig. [21 is a particularly simple, one parameter species and there is no reason 
to expect that it is already optimal. Rather, paths realized in experiments 
will most likely be somewhat different, i.e. some (smooth) deviations of it. 
(3) It is a well known fact, that due to the local approximations in the ex- 
change correlation functionals (XC), DFT calculations do not usually give 
precise quantitative results for conductances [32] or molecular polarizabili- 
ties [23]. For this reason, our theoretical results have an intrinsic quantitative 
uncertainty. 

We believe, that in view of points (1-3) it is not surprising that our 
estimate for the switching value Kn~2— 4 V does not very precisely reproduce 
the experimental one ~ 1 V [T^ |2I] . 

8 Predictions 

We close with offering two qualitative predictions that follow from our theo- 
retical analysis: (pi) The "easy axis" rotational degree of freedom of BPDN 
around the carbon triple bonds. Fig. [H is crucial for 6'-rotations. If easy rota- 
tion is blocked e. g. by molecular design, then the rotational barriers become 
much higher, switching is more difficult and hysteresis should be suppressed. 
(p2) The switching voltage. Von can be modified by dressing the BP-unit, 
Fig. dl with other redox groups, e.g. carboxy-groups, COOH, instead of the 
nitro-groups, NO2. We expect, that the switching voltage decreases when the 
intrinsic dipole moment of the BP-complex becomes larger. (Here, we have 
assumed that the change in molecular conformation, which accompanies the 
exchange of the redox groups, is less important.) 

Summarizing, we have identified a two-axis rotation that takes a Tour- 
type, BPDN-DT, molecular junction from its ground-state conformation 
("off") into a metastable one ("on"). The energy barrier for this process 
is larger than temperature, but can be controlled by applying a bias volt- 
age. At sufficiently large biases, the molecule undergoes the transition from 
"off" -state to "on" and reverse, so that it can realize a molecular switch. Our 
theoretical analysis may explain the pronounced switching behavior observed 
in recent experiments [T9l [20l [2T] . 
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BPDN-DT 

O2N 




NO2 

BP-DT 




Figure 1: Upper panel: The bipyridyl-dinitro oligophenylene-ethynylene 
dithiol (BPDN-DT) molecule that shows switching in single molecule trans- 
port experiments [231 1201 121] • The protection groups, Ac, are released when 
the molecule forms a chemical bond with the Au-surface. Lower panel: a 
reference molecule that does not switch [20l |2T]. 
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Figure 2: Extended molecule used for the DFT calculations. The NO2 groups 
(red/blue) introduce a dipole moment perpendicular to the wire direction. 
In addition, their steric interaction drives the wire out of the planar confor- 
mation. 
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Figure 3: Energies of the highest occupied (HOMO, sohd, red) and the low- 
est unoccupied (LUMO, dashed, green) molecular levels, for the molecular 
wire depicted in Fig. [1] (with the replacement Ac— s>H) . Also neighboring or- 
bitals, HOMO-1 (solid, brown) and LUMO+1 (dashed, blue), are shown. For 
the molecule with excess charge, Q= ± 1, orbitals have been obtained with 
( "screened" ) and without ( "unscreened" ) image charges; this models screen- 
ing in the vicinity of metallic (Au) electrodes. Fermi energy of Au fcc-clusters 
are i?Fermi~ — 5.05eV, with BP-functional (cf. i?Fermi= — 5.65, —5.55, —5.42 eV 
for Au(lll), Au(lOO), Au(llO) with local density approximation (LDA) [HT] . 
respectively). The plot on the left (a) shows results obtained within the 
BP86 [28l [29] exchange correlation functonal, while the plot of the right (b) 
shows corresponding results for the B3LYP [30] functional. No qualitative 
difference in level positions is observed indicating that self-interaction errors 
do not interfere. 
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Figure 4: Schematics indicating the two axis rotation, which connects the two 
bistable states reahzing the molecular switch. Blue arrow indicates the NO2 
group, grey lines the carbon backbone, cf. Fig. [2l First step is rotation of 
the whole molecule by 6 about the red axis connecting two S-atoms. Second 
step is (back-)rotation of the connecting phenyl ring only about the green axis 
with angle —9. Back rotation partly eliminates the energy cost associated 
with contact modification. Grey arrow indicates the direction of an external 
electric field. 
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Figure 5: Dependency of the total ground state (DFT) energy, AEtot{0), of 
the extended molecule (Aui4-clusters) on the rotation angle 6 and an inter- 
electrode electric potential drop, V, with homogeneous gradient E = — VV^, 
see Fig. HI (The simulation angle is $=70°; full geometry relaxation has been 
done for 1/ = V at ^ = 0.) 
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Figure 6: Current voltage trace for the BPDN-DT molecule at $^70° for 
the ground state configuration, GC (o), as well as for the meta-stable one, 
MC (□). Inset upper left: blow up of main plot near small negative Vbias- 
Inset lower right: Behavior of IV curves under the transformation Vbias 
-Vbias, / ^ dashed line (left "on", right "off"). 
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